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Abstract 

Fractional derivative can be denned as a fractional power of derivative. The 
commutator (i/h)[H, . ], which is used in the Heisenberg equation, is a derivation 
on a set of observables. A derivation is a map that satisfies the Leibnitz rule. In 
this paper, we consider a fractional derivative on a set of quantum observables as a 
fractional power of the commutator (i/h)[H, .]. As a result, we obtain a fractional 
generalization of the Heisenberg equation. The fractional Heisenberg equation is 
exactly solved for the Hamiltonians of free particle and harmonic oscillator. The 
suggested Heisenberg equation generalize a notion of quantum Hamiltonian systems 
to describe quantum dissipative processes. 

PACS: 03.65.-w; 03.65. Ca; 45.10.Hj; 03.65.Db 

Keywords: Heisenberg equation, fractional derivative, fractional equation 



1 



1 Introduction 



The fractional calculus has a long history from 1695, when the derivative of order 1/2 has 
been described by Leibniz [TJ E] . The theory of derivatives and integrals of non- integer 
order goes back to Leibniz, Liouville, Grunwald, Letnikov and Riemann. There are many 
books about fractional calculus and fractional differential equations [TJ El El H] . Derivatives 
of fractional order, and fractional differential equations have found many applications in 
recent studies in physics (see for example [HEJITIIE] and references therein). 

In the quantum kinematics, the observables are given by self-adjoint operators. The 
dynamical description of system is given by a superoperator, which is a rule that assigns 
to each operator exactly one operator. The natural description of the motion is in terms of 
the infinitesimal change of the system. The equation of motion for quantum observable is 
called the Heisenberg equation. For Hamiltonian systems, the infinitesimal superoperator 
of the Heisenberg equation is defined by some form of derivation. 

Fractional derivative can be defined as a fractional power of derivative (see for exam- 
ple [H]). It is known that the infinitesimal generator (i/h)[H, . ], which is used in the 
Heisenberg equation, is a derivation of observables. A derivation is a linear map D, which 
satisfies the Leibnitz rule D(AB) = (DA)B + A(DB) for all operators A and B. In this 
paper, we consider a fractional derivative on a set of observables as a fractional power of 
derivative. As a result, we obtain a fractional generalization of the Heisenberg equation. 
It allows us to generalize a notion of quantum Hamiltonian systems. Note that fractional 
generalization of classical Hamiltonian systems has been suggested in [TU] . The suggested 
fractional Heisenberg equation is exactly solved for the Hamiltonians of free particle and 
harmonic oscillator. A quantum system that is presented by fractional Heisenberg equa- 
tion can be considered as a dissipative system [IT] . Fractional derivatives can be used as 
a possible approach to describe an interaction between the system and an environment. 
Note that it is possible to consider quantum dynamics with low-level fractionality by some 
generalization of method suggested in [12] (see also [13]). 

In Section 2, the fractional power of derivative and the fractional Heisenberg (FH) 
equation are suggested. In Section 3, the properties of time evolution, which is described 
by the fractional equation, are considered. In Section 4, the FH equation for free particle 
is solved. In Section 5. the solution of FH equation with harmonic oscillator Hamiltonian 
is derived. 
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2 Fractional derivative and Heisenberg equation 



Let us consider a set of quantum observables in the Heisenberg picture. A superoperator £ 
is a rule that assigns to each operator A exactly one operator £(A). (About superoperator 
formalism see for example [T4"l [151 [16] and references therein.) For Hamiltonian H, let 
L^j be the superoperator given by 

L~ H A = \[HA- AE). 
in 

The operator differential equation 

iA = - L « A ' w 

is called the Heisenberg equation for Hamiltonian systems. The time evolution of a 
Hamiltonian system is induced by the Hamiltonian H. 

It is interesting to obtain a fractional generalization of equation (CQ). We will consider 
here concept of fractional power for Ljj. If is a closed linear superoperator with an 
everywhere dense domain D(L]j), having a resolvent R(z,Ljj) = (zLj — Ljj)^ 1 on the 
negative half-axis, then there exists [19], [20j [21] the superoperator 

- (L H ) a = r dz z«-'R(-z, L H ) L H (2) 

7T JO 

defined on D(L]j) for < a < 1. The superoperator is a fractional power of the 

Lie left superoperator. 

As a result, we obtain the equation 

jA t = -{L H YA U (3) 

where t and H/h are dimensionless variables. This is the fractional Heisenberg equation. 
Note that this equation cannot be presented in the form 

j t A = -L Hnew A t = ^-[H new ,A t ] (4) 

with some operator H new . Therefore, quantum systems described by ([3]) are not Hamil- 
tonian systems. The systems will be called the fractional Hamiltonian systems (FHS). 
Usual Hamiltonian systems can be considered as a special case of FHS. 

If we consider the Cauchy problem for equation ([T]) in which the initial condition is 
given at the time t = by A , then its solution can be written in the form 
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The one-parameter superoperators $(, t > have the properties 



$i$ s = $t +s , (t,s>0), % = L r , 

where Lj is unit superoperator [LiA = A). Then the set of $t, t > 0, is called the semi- 
group. Then the superoperator is called the generating superoperator, or infinitesimal 
generator, of the semi-group {$i,t > 0}. 

Let us consider the Cauchy problem for fractional Heisenberg equation ([3]) in which 
the initial condition is given by Aq. Then its solution can be presented in the form 

A t (a) = $| a) A , 

where the superoperators t > 0, form a semi-group which will be called the frac- 
tional semi-group. The superoperator {L^) 01 is infinitesimal generator of the semi-group 
{$| Q) ,t>0}. 

3 Properties of fractional time evolution 

Let us consider some properties of time evolution described by a fractional semi-group 
{$| Q) ,t>0}. 

(1) The superoperators <£>^ can be constructed in terms of the superoperators $t by 
the Bochner-Phillips formula [TT1 [TS1 12U] 

= rdsf a (t,s)$ s , (t>0), (5) 
Jo 

where f a (t,s) is defined by 

fa(t, s) = — / dz exp(s2; - tz a ), (6) 

where a, t > 0, s > 0, and < a < 1. The branch of z a is so taken that Re(z a ) > for 
Re(z) > 0. This branch is a one-valued function in the z-plane cut along the negative 
real axis. The convergence of this integral is obviously in virtue of the convergence 
factor exp(— tz a ). By denoting the path of integration in ([6]) to the union of two paths 
r exp(— i8), and r exp(+i6 l ), where r G (0, oo), and tt/2 < 9 < n, we obtain 

1 r°° 

f a (t,s) = — dr exp(sr cos 9 — tr a cos(a#))- 
7T Jo 

• sin(sr sin 9 — tr a sin(a^) + 9). (7) 
If we have a solution A t of the Heisenberg equation ([ID, then formula © gives the solution 

POO 

A t (a) = / dsf a {t,s)A s , (t>0) (8) 
Jo 
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of fractional equation (Ej). 

(2) In quantum theory, the most important is the class of real superoperators. Let 
A* be an adjoint of A. A quantum observable is a self-adjoint operator. If $ t is a real 
superoperator and A is a self-adjoint operator A* = A, then the operator A t = & t A is 
self-adjoint, i.e., 

($ t A)* = $ t A 

A superoperator as a map from a set of observables into itself should be real. All possi- 
ble dynamics, i.e., temporal evolutions of quantum observables, should be described by 
real superoperators. Therefore the following statement is very important. If $ t is a real 
superoperator, then <3>j Q ' ) is real. It can be proved by using the Bochner-Phillips formula 
and equation (J7J). 

(3) An adjoint superoperator of $ t is a superoperator 3> 4 , such that 

Tr[($ t A)*B] = Tr[A*$ t (B)]. 

Let us give the basic statement regarding the adjoint superoperator. If $ t is an adjoint 
superoperator of <E>t, then the superoperator 



$i a) = rdsf a (t, s )$ s , (t>o), 

Jo 



is an adjoint superoperator of ■ We prove this statement by using the Bochner-Phillips 

formula: ^ 

Tr[(^ a) A)*B] = f dsf a (t,s)Tr[(® s A)*B] = 
Jo 



rdsf a (t,s)Tr[A*$ s (B)]=Tr[A*$i a) (B)\. 
Jo 



Let {$t,t > 0} be a semi-group, such that the density matrix operator p t = $ t p is 
described by the von Neumann equation 

Then the semi-group {^[ a \t > 0} describes the evolution of the density operator 

p t (a) = $| Q) po 

by the fractional equation 

j t Pt = -(-L H ) a p t . 
This is the fractional von Neumann equation. 
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(4) Let $t, t > 0, be a positive one-parameter superoperator, i.e., $fv4 > for A > 0. 
Using the Bochner-Phillips formula and the property f a (t,s) > 0, (s > 0), it is easy to 
prove that the superoperators $^ are also positive, i.e., $\ a A > for A > 0. 

(5) It is known that $ t is a real superoperator if $ t is real. Analogously, if is a 
real superoperator, then $1°^ is real. 



4 Fractional free particle 

Let us consider the Hamiltonian H = P 2 /2m, where P is dimensionless variable, and 
m _1 has the action dimension. Then the Heisenberg equation ([I]) describes a free one- 
dimensional particle. For A = Q, and A = P, equation ([I]) gives 

d „ 1 c2 

-Q t = -P*, -P = 0. 

etc m at 

The well-known solutions of these equations are 

Qt = Qo + -Po, Pt = Po- (9) 
m 

Using these solutions and the Bochner-Phillips formula, we will obtain the solutions of 
fractional Heisenberg equations 



|ft = ~(^m, |* = o. do) 



Note that (L p2 ) a ^ L p2a - The solutions of fractional equations (fTUl) have the forms 

Q t (a) = $S a) Q = [°° dsf a (t, s)Q 8 , P t (a) = P , 
Jo 

where Q s is given by ©. Then 

= Qo + — a (t)Po, ^ = Po, 

m 

where _ 



9a(t) = I dsf a (t, S) S. 

If a = 1/2, then we have 



o 



9i/2 (t) = —= / ds—=e > is = -. 
2a/7t Jo y/s 2 
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Then 2 

Qt = Qo- ^-Po, Pt = Po- (11) 

2m 

These equations describe a fractional free motion for a = 1/2. For the operators f fTTl) . 
the average values and dispersions have the form 

t 2 

< Q t >= x - —Po, < Pt >= Po, 
2m 

and 

Here we use the coordinate representation and the pure state 

V(x) =< x\V >= (fe v ^)- 1 / 2 exp(- (x ~* o) * + ip z). (12) 

The average value and dispersion are defined by the well-known equations 

< A t >= Tr[|* >< *|A t ] =< V\A t \y >, 
Dj^it) =< A 2 > - < A* > 2 =< ^K 2 |* > - < *|A t |* > 2 . 

5 Fractional Heisenberg equation for harmonic oscil- 
lator 

Let us consider the Hamiltonian 

1 9 muj 2 ~ 
H = P 2 + Q 2 , 13 

where £ and P are dimensionless variables. Then equation describes a harmonic 
oscillator. For A = Q, and A = P, equation flTJ gives 

^Q* = -P;, ^P t = -m^ 2 Q,. 
at m at 

The well-known solutions of these equations are 

Qt = Qo cos(urb) H P sin(co>t), 

Pt = Pocos(cijt) — mujQo sin(a>t). (14) 
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Using these solutions and the Bochner-Phillips formula, we will obtain the solutions of 
fractional Heisenberg equations 



j t Q t = -(L H ) a Q t , jP t = -{L H YP U 



(15) 



where H is defined by (1131) . The solutions of fractional equations (|T5|) have the forms 

Q t (a) = $| a) Q = r > dsf a (t,s)Q s , 
Jo 



P t (a) = $| a) P 



dsf a (t,s)P a . 



Substitution of (TH|) into (|T6|) gives 

Qt — QoC a (t) 



l 



PoS a (t), 



P t = P C a (t)-muQ S a (t), 



(16) 

(17) 
(18) 



where 



C a (t) 
S a {t) 



ds fait, s) cos(o;s), 
ds f a {t, s) sin(o;s). 



Equations (TTTT) and (|T8|) describe solutions of the fractional Heisenberg equations (fl5T) for 
quantum harmonic oscillator. 
If a = 1/2, then 



C V 2{t) 
Sx/2(t) 



t 



2^ Jo 
t 



ds 

i 

ds 



cos ws 



s 3/2 

x sin(o;s) 



-t 2 /4s 



-t 2 /4s 



,3/2 



These functions can be presented through the Macdonald function (see [22], Sec. 2.5.37.1.) 
such that 



C 1/2 (t) = (^) 1/4 [e + - /8 ^- 1/2 (2e + -/ 4 V /^) +e-^_ 1/2 (2 C -W^) 



f 2 - 



/u^_ 2 - 

T 



= <17) 1/4 [ e+W ^-V 2 ( 2e+W Vx) - e ^ /8 ^i/ 2 (2e- /4 ^)], (19) 

where ui > 0, and lf a (z) is the Macdonald function [HE], which is also called the modified 
Bessel function of the third kind. 



tut 2 
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Using (TI2I) . we get the average values 



< Q t >= x C a (t) H PoS a (t), 

moo 

< P t >= p C a (t) - mwx S a (t), 

and the dispersions 

Dp(t) = ^ 2 Cl(t) + b ^^Sl(t), 
2 2o z m z uj z 

As a result, the fractional harmonic oscillator is a simple dissipative systems. The disper- 
sion of the wave packet is defined by these equations. The solutions are characterized by 
the fractional damping of the fractional harmonic oscillator. The dumping is described 
by the modified Bessel function ( TT9|) of the third kind. An important property of the 
evolution described by the fractional equations are that the solutions have power-like 
tails. 



6 Conclusion 

In this paper, we consider derivatives of noninteger order as fractional powers of derivative. 
We derive a fractional generalization of the Heisenberg equation and a generalization of 
quantum Hamiltonian system to describe open quantum systems. A quantum system that 
is presented by fractional Heisenberg equation can be considered as a dissipative system. 
Fractional derivatives can be used as an approach to describe an interaction between the 
quantum system and an environment. This interpretation caused by following reasons. 
Using the properties 

POO 

/ f a (t,s) = l, f a (t,s)>0 (for all s > 0), 
Jo 

we can assume that f a (t,s) is a density of probability distribution. Then the Bochner- 
Phillips formula © can be considered as a smoothing of Hamiltonian evolution $ t with 
respect to time s > 0. This smoothing can be considered as an influence of the environ- 
ment on the system. As a result, the parameter alpha can be used as a simple model 
of the interaction between the system and the environment. Note that the quantum 
Markovian equations, which is also called the Lindblad equations [23J , can be generalized 
by suggested approach to describe completely positive evolution of dissipative and open 
quantum systems. 
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